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Abstract 

We demonstrate that in a coupled two-qubit system any single-qubit gate can be decomposed into two conditional two-qubit 
gates and that any conditional two-qubit gate can be implemented by a manipulation analogous to that used for a controlled 
two-qubit gate. Based on this we present a unified approach to implement universal single-qubit and two-qubit gates in a 
coupled two-qubit system with fixed always-on coupling. This approach requires neither supplementary circuit or additional 
physical qubits to control the coupling nor extra hardware to adjust the energy level structure. The feasibility of this approach 
is demonstrated by numerical simulation of single-qubit gates and creation of two-qubit Bell states in rf-driven inductively 
coupled two SQUID flux qubits with realistic device parameters and constant always-on coupling. 

PACS numbers: 03.67.Lx, 85.25.Dq, 03.67.Mn 
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During the past decade, a variety of physical qubits 
have been explored for possible implementation of quan- 
tum gates. Of those solid-state qubits based on su- 
perconducting devices have attracted much attention 
because of their advantages of large-scale integration 
and easy connection to conventional electronic circuits 
Hiiaiillli M MMMMM -Super- 
conducting single-qubit gates H, ||J la 0, la M, ilii LLil an d 
two-qubit gates 0, 0, 0, llfij have been demonstrated 
recently. 

However, building a practical quantum computer re- 
quires to operate a large number of multi-qubit gates si- 
multaneously in a coupled multi-qubit system. It has 
been demonstrated that any type of multi-qubit gate 
can be decomposed into a set of universal single-qubit 
gates and a two-qubit gate, such as the controlled-NOT 
(CNOT) gate Hi E3- Thus it is imperative to imple- 
ment the universal single-qubit and two-qubit gates in a 
multi-qubit system with the minimum resource and max- 
imum efficiency |l8j . 

Implementing universal single-qubit gates and two- 
qubit gates in coupled multi-qubit systems can be 
achieved by turning off and on the coupling between 
qubits ^1, ElUcJ- In these schemes, supplementary cir- 
cuits were required to control inter-qubit coupling. How- 
ever, rapid switching of the coupling results in two serious 
problems. The first one is gate errors caused by popu- 
lation propagation between qubits. Because the compu- 
tational states of the single-qubit gates are not a subset 
of the eigenstates of the two-qubit gates the populations 
of the computational states propagate from one qubit to 
another when the coupling is changed, resulting in addi- 
tional gate errors. The second one is additional decoher- 
ence introduced by the supplementary circuits 0, H^] • 
This is one of the biggest obstacles for quantum comput- 
ing with solid-state qubits, particularly in coupled multi- 
qubit systems Q. In addition, the use of supplementary 
circuits also significantly increase the complexity of fab- 



rication and manipulation of the coupled qubits. 

To circumvent these problems, a couple of alternative 
schemes, such as those with untunable coupling |2l| and 
always-on interaction (2^, have been proposed. In the 
first scheme, each logic qubit is encoded by extra physical 
qubits and coupling between the encoded qubits is con- 
stant but can be turned off and on effectively by putting 
the qubits in and driving them out of the interaction free 
subspace. In the second scheme, the coupling is always 
on but the transition energies of the qubits are tuned in- 
dividually or collectively. These schemes can overcome 
the problem of undesired population propagation but 
still suffer from those caused by the supplementary cir- 
cuits needed to move the encoded qubits and tune the 
transition energies. Moreover, the use of encoded qubits 
also requires a significant number of additional physical 
qubits. 

In this Letter, we present a unified approach to imple- 
ment universal single-qubit and two-qubit gates in a cou- 
pled two-qubit system with fixed always-on coupling. In 
this approach, each single-qubit gate is realized via two 
conditional two-qubit gates and each conditional two- 
qubit gate is implemented with a manipulation analogous 
to that used for a controlled two-qubit gate (e.g., CNOT) 
in the same subspace of the coupled two-qubit system 
without additional circuits or physical qubits. Since the 
computational states of the single-qubit gates are a sub- 
set of those of the two-qubit gates the gate errors due to 
population propagation are completely eliminated. The 
effectiveness of the approach is demonstrated by numer- 
ically simulating the single-qubit gates and creating the 
Bell states in a unit of inductively coupled two supercon- 
ducting quantum interference device (SQUID) flux qubits 
with realistic device parameters and constant always on 
coupling. 

Consider a basic unit consisting of two coupled qubits 
which we call a control qubit (Cq) and a target qubit (Tq) 
for convenience. An eigenstate of the coupled qubits is 
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denoted by |n) = \ij) , which can be well approximated by 
the product of an eigenstate of Cq, and that of Tq, 
\n) = \ij) = \i) for weak inter-qubit coupling. The 
computational states of the coupled qubits are |1) = 1 00} , 
1 2) = |01), 1 3) = 1 10), and |4) = |11), which correspond 
to the eigenstates for i, j = 0, 1 and n = 1, 2, 3, 4, 
respectively. In general, the result of an operation on 
Cq depends on the state of Tq since Cq is coupled to 
and hence influenced by Tq. Suppose U is a unitary 
operator acting on a single qubit in the four-dimensional 
(4D) Hilbert space spanned by \n). In order to perform 
U on Cq one wants the state of Cq evolves according to U 
independent of the state of Tq which is left unchanged. 
This operation is denoted by = ® , where, 
the 4x4 unitary matrix and the 2x2 unitary matrix 



U. 



q are the representations of U in the Hilbert space 
of the coupled qubits and in the subspace of Cq while 
the 2x2 unitary matrix lQ' is the representation of the 
identity operation in the subspace of Tq. If the matrix 
elements of are u%j = 1, 2) the operation can be 
decomposed into two operations as 
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where, the 4x4 unitary matrices U^ 1 and U[ are also 
in the 4D Hilbert space of the coupled qubits. Since 
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only involves the states |00) and 1 10} and U^' only 
involves |01) and |11) they represent operations on Cq 
when Tq is in the state |0) and |1), respectively. If the 

(2) (2) 

state of Tq is |0) (|1)), Uq (E/f ) implements the desired 

(2) (2) (2) 

operation U while U{ (Uq ) does nothing. Thus Uq 

(2) 

and U{ represent two conditional two-qubit gates of the 
coupled qubits. Eq.QJ shows that for the coupled qubits 
any single-qubit gate can be realized via two conditional 
two-qubit gates. This property is universal and is shown 
schematically in Fig. 1, where the single-qubit operation 
C/( 2 ) and its equivalent operation via two conditional two- 

(2) (2) 

qubit gates Uk and U{ are illustrated. The conditional 
gates are denoted by the open and closed circles for the 
state of Tq being |0) and |1), respectively [lij . 

The most obvious advantage of implementing single- 
qubit gates this way is that the conditional two-qubit 
gates can be realized using essentially the same set of 
operations used for controlled two-qubit gates. To show 
this, let us consider a rf-driven coupled two-qubit system 
with sufficiently different energy- level spacings AE13, 
AE24, Ai?i2, and AE34, where AE nn ' is the level spacing 
between the states \n) and \n'). To implement, for exam- 
ple, a single-qubit NOT gate, we apply two n pulses to 
the Cq: the first one is resonant with AE13 and the second 
with Ai?24- Because both pulses are largely detuned from 
AE12 and Ai?34, they do not induce unintended popula- 
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FIG. 1: Single-qubit gate and its equivalent operation 

(2) (2) 

by two conditional two-qubit gates Uq and U{ in a coupled 
two-qubit system. 



tion transfer when the fields are sufficiently weak. If the 
state of Tq is |0) (|1)), the first (second) pulse accom- 
plishes the desired transformation. Therefore, when the 
two microwave pulses are applied to Cq, either sequen- 
tially or simultaneously, the NOT gate is accomplished 
and the gate time is essentially the same as that of a 
stand-alone Cq if the both pulses are applied simultane- 
ously. 

To implement a two-qubit gate such as CNOT in cou- 
pled qubits, we apply a ir pulse resonant with A-E34 to Tq 
|24j . In this case, the state of Tq flips if and only if the 
state of Cq is |1). Since the energy level structure for the 
conditional two-qubit gates and the CNOT gate could 
be the same the universal single-qubit gates and CNOT 
gate can be implemented in the same fashion without 
adjusting inter-qubit coupling as long as Cq and Tq can 
be addressed individually by microwave pulses, which is 
rather straightforward to realize with solid-state qubits 
in general and with rf SQUID flux qubits in particular. 

For concreteness, we demonstrate how to realize the 
single-qubit and two-qubit gates in rf-driven two SQUID 
flux qubits with constant always-on coupling. The cou- 
pled flux qubits comprise two SQUIDs coupled induc- 
tively through their mutual inductance M. Each SQUID 
consists of a superconducting loop of inductance L in- 
terrupted by a Josephson tunnel junction characterized 
by its critical current I c and shunt capacitance C |2fil |. 
A flux-biased SQUID with total magnetic flux $ en- 
closed in the loop is analogous to a "flux" particle of 
mass m = C'&q moving in a one-dimensional poten- 
tial, where $0 = h/2e is the flux quantum. For sim- 
plicity, we assume that the two SQUIDs are identical. 
In this case, Cj = C, = L, and I C i = I c for 
i = 1.2. The Hamiltonian of the coupled qubits is 
013 H{ Xl ,x 2 ) = H ( Xl ) +H {x 2 ) +H 12 (x u x 2 ), 
where Hq (x^) is Hamiltonian of the zth single qubit given 
by H (xi) = Pt/2m+mu1 c (x t - x ei ) /2—Ej cos (2tt Xi ) 
and H12 is the interaction between the SQUIDs given 
by i?i2 = mu)\ c K{xx - x e \) (x 2 - x e2 ) /2. Here, x t = 
$i/&o is the canonical coordinate of the ith "flux" parti- 
cle and pi = —ihd/dxi is the canonical momentum con- 
jugate to Xi, x e i = <& e i/<f>o is the normalized external 
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flux bias of the ith qubit, Ej — mw| /C ,/3i/47r 2 is the 
Josephson coupling energy, /3l — 2irLI c /$ is the poten- 
tial shape parameter, lolc = 1/VLC is the characteristic 
frequency of the SQUID, and K = 2M/L is the coupling 
strength. The coupled SQUID qubits are a multi-level 
system. The eigenenergies E n and eigenstates \n) are ob- 
tained by numerically solving the eigenvalue equation of 
H(x\,X2) using the two-dimensional Fourier-grid Hamil- 
tonian method When x e \ and x e 2 ~ 0.5 the coupled 
SQUID qubits have four wells in the potential energy sur- 
face [24$ . The four computational states are chosen to be 
the lowest eigenstate of each well. 

To realize single-qubit and two-qubit gates in 
the coupled SQUID qubits, we apply microwave 
pulses xc and xt to Cq and Tq, respectively. 
The interaction of the coupled qubits and pulses is 
V (x\,X2 1 t) — d\ (xi — x e \) + di (x 2 — x e2 ) + du with 



do 



(xc + kxt/2), d 2 = 



J LC 
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(x c + x T + kxcxt) /2. 



kxc/2), 
The 



d\ = mull 
and d\2 

time-dependent wave functions of the coupled 
SQUID qubits, ip(x\, X2, i), are governed by the 
time-dependent Schrodinger equation ihdip/dt — 
[H(xi, X2) + V(xx, X2, t)] ip. To solve this equation we ex- 
pand ip in terms of the first 20 eigenstates of the coupled 
qubits: ip — ^2 n c n (t) \n). The expansion coefficients 
c n (r) are calculated by numerically solving the matrix 
equation i<9c„ (r) jdr = J2 n > H nn> ( T ) c «' ( T ) usin S tne 
split-operator method |28| . where r = u>Lct is the dimcn- 
sionless time and H^ n , = [E n S nn ' + (n \V\ n')} /fkVLC is 
the reduced Hamiltonian matrix element. The probabil- 
ity of being in the state \n) is thus |c„(t)| 2 . For single- 
qubit gates two resonant pulses, xci — xcw cos (u>cit) 
with ujcx — AEi3/h and xc2 = XC20 cos (u>c2t) with 
UJC2 — A£ , 24/S, are applied simultaneously to the Cq 
so that xc = xci + xc2- For controlled two-qubit 
gates such as the CNOT gate one resonant pulse 
xt = xto cos (ojxt) with lot = AE^/h is applied to 
the Tq. To minimize the possible intrinsic gate errors 
caused by leakage to unintended states, we select the 
parameters of the coupled SQUID qubits using the 
independent transition approach [2^|. For SQUIDs with 
L = 100 pH, C = 40 fF, and f3 L = 1.2 one set of the 
better working parameters is x e \ — 0.499, x e 2 = 0.4998, 
and k = 5 x 10~ 4 , which will be used in the calculation. 

Most of the single-qubit gates, such as the NOT gate 
and Hadamard gate, can be described by rotations of 
the state vector on the Bloch sphere [la]. Thus we 
demonstrate how to realize an arbitrary single-qubit ro- 
tation of an angle 9 about an axis perpendicular to the 
z axis, RS 2 '(6), in the coupled SQUID qubits. Based on 
Eq.lJTJ, the rotation R^ 2 \9) on Cq is accomplished via 

(2) (2) 

two conditional two-qubit rotations Rq (9) and R\ ' (9) 
as shown in FIG. 2(a). They are realized by apply- 
ing two microwave pulses xci and xc2 to Cq simulta- 
neously. The amplitudes and frequencies of xci and xc 2 
are x C io = 5 x 10" 5 , w C i = 0.239w ic = 2ir x 19.0 GHz, 
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FIG. 2: Arbitrary single-qubit rotation about an axis perpen- 
dicular to z axis in the coupled SQUID flux qubits. (a) Quan- 
tum circuit: the state \ij) evolves to the state RP'(6) \ij) = 
^R( 2 \9) \j) after the two conditional two-qubit rotations. 

(b) The rotation angle 9 vs. pulse width, (c) Populations of 
the states |00) and 1 10) vs. pulse width. 



XC20 = 5.14 x 10~ 5 , and lj C 2 = 0.259w LC = 2ir x 20.6 
GHz. In Fig. 2 (b) and (c), we plot the rotation angle 
9 and the populations of the states 1 00) and 1 10) (pop- 
ulations of the other states remain essentially at zero) 
as a function of pulse width when the initial state of 
the coupled qubits is |00). It is shown in Fig. 2 (b) 
that the rotation angle 9 is essentially a linear function 
of pulse width. This indicates that the state of Cq un- 
dergoes Rabi oscillations for which the phase angle fir 
is a linear function of pulse width, where f2 is the Rabi 
frequency. It is also shown in Fig. 2 (c) that from the 
initial state 1 00) the coupled qubits evolve into the state 
( 1 00) + 1 10)) /v2 after the tt/2 pulses and into the state 
1 10) after the it pulses. We have also computed the rota- 
tion angles and populations for the coupled qubits with 
the initial states 1 01) , 1 10) , and |11) using the same pulse 
sequence. In each case, the state is transformed from 
\ij) to RW(6) \ij) = [RW(d) \j) for i,j = 0,1 with 
6 = Ht as expected. 

Entanglement is one of the most profound charac- 
teristics of quantum systems which plays a crucial role 
in quantum information processing and communication 
[lfj| . The maximally entangled two-qubit states are re- 
ferred to as the Bell states. To create the Bell states 
from a state \ij), a Hadamard gate on Cq, H^ 2 \ which 

is decomposed into two conditional two-qubit Hadamard 

(2) (2) 

gates Hq and H\ , and a CNOT gate are commonly 
used [see FIG. 3(a)]. The two conditional two-qubit 
Hadamard gates are implemented by applying two tt/2 
pulses xci and xc2 to Cq and the following CNOT gate 
is implemented by applying a it pulse xt to Tq, as shown 
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FIG. 3: Creation of the Bell states in the coupled SQUID flux 
qubits. (a) Quantum circuit: the state \ij) evolves into an en- 
tangled state \ Bij) after two conditional two-qubit Hadamard 
gates and one CNOT gate, (b) Microwave pulses: two n/2 
pulses xci and xq2 are applied to the Cq first and then a n 
pulse xt is applied to the Tq. (c) Population evolution: the 
coupled qubits evolve into a product state ( 1 00) + 1 10} ) /a/2 
from the initial state 1 00) after the first two n/2 pulses and 
then into a Bell state ( 1 00} + |11)) /v2 from the product state 
after the second it pulse. 



in FIG. 3(b). The amplitudes and frequencies of xci and 
xc2 are the same as those used in FIG. 2 and those of 
x T are x T q — 5 x 10~ 5 and lu t = 0.0592w ic — 2n x 4.7 
GHz. In FIG. 3(c), we plot the population evolution of 
the computational states when the initial state is 1 00) . 
It is shown clearly that the coupled qubits evolve first 
into a product state ( 1 00) + 1 10) ) /V2 from the initial 
state 1 00) after the ir/2 pulses and then into a Bell state 
( 1 00) + 1 1 1) ) /\/2 after the subsequent tt pulse. We have 
also calculated the population evolution for the coupled 
qubits being initially in 1 01) , |10), and |11), respectively. 
The final state in each case is also one of the expected 
Bell states. 

In summary, we showed that in a coupled two-qubit 
system any single-qubit gate can be realized via two condi- 
tional two-qubit gates and that any conditional two-qubit 
gate can be implemented with a manipulation analogous 
to that used for a controlled two-qubit gate. Based on 
this universal property of single-qubit gates we present 
a general approach to implement the universal single- 
qubit and two-qubit gates in the same coupled two-qubit 
system with fixed always-on coupling. This approach is 
demonstrated by using a unit of two SQUID flux qubits 
with realistic device parameters and constant always-on 
coupling. Compared to other methods our approach has 
the following characteristics and advantages: (1) The ap- 
proach is universal as long as each qubit can be locally 



addressed; (2) No additional decoherence from the hard- 
ware added to control inter-qubit coupling; (3) Gate error 
induced by the population propagation from one qubit 
to another is completely eliminated; (4) The architec- 
ture for both hardware (circuits) and software (pulse se- 
quence) is much simplified. This approach can be readily 
extended to multi-qubit systems and other types of solid- 
state qubit systems in which each qubit can be individu- 
ally addressed. Therefore, it is very promising for realiz- 
ing universal gates with minimum resource and complex- 
ity and maximum efficiency. 

This work is supported in part by the NSF (DMR- 
0325551) and by AFOSR, NSA, and ARDA through 
DURINT grant (F49620-0 1-1-0439). 
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